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Abstract 

In this paper we study biased random K-SAT problems in which each logical variable is negated 
with probability p. This generalization provides us a crossover from easy to hard problems and 
would help us in a better understanding of the typical complexity of random K-SAT problems. The 
exact solution of 1-SAT case is given. The critical point of K-SAT problems and results of replica 
method are derived in the replica symmetry framework. It is found that in this approximation 
a c oc for p —* 0. Solving numerically the survey propagation equations for K = 3 we find 

that for p < p* ~ 0.17 there is no replica symmetry breaking and still the SAT-UNSAT transition 
is discontinuous. 
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I. INTRODUCTION 



Optimization problems are subject of recent studies in the context of complex systems 



flung 



Random K-SAT problems are well known examples of these problems which 
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have their origin in computer science and complexity theory [5|, |6|, 13, |8|, |9|, |10|, |11 1 . Finding 
the configuration of N logical variables which satisfies a formula of M clauses is a hard 
problem and indeed lies in the class of NP-complete problems for K > 3^^. 
From a physical point of view the interesting feature of random K-SAT problems is the 
presence of phase transitions in the thermodynamic limit where N and M approach 
infinity and a := M/N remains finite |2|. Here the transition is between SAT and UNSAT 
phases where a typical instance of the formula is satisfied or unsatisfied respectively with 
probability 1. It is around the phase transition that the time needed to find the solution of 
a typical instance grows exponentially with the size of the problem N. 
In this paper we study a generalized version of the random K-SAT problems where each 
logical variable is negated with probability p rather than 1/2 as in the original random 
K-SAT. The aim of this generalization is to go continuously from easy instances of the 
problem to hard ones. Clearly for p = we have an easy random K-SAT for all values 
of K. On the other hand as the studies indicate, the problem is hard for p = 1/2 and 
K > 3. Thus one expects a crossover from easy to hard region by increasing p from zero. 
Is it possible to define a point beyond which one can say that the problem becomes hard? 
How the problem approaches the hard regime? What are the universal features of this 
crossover? These are some questions which can provide us a deeper understanding of the 
typical complexity of random K-SAT problems. 

A similar problem to one that we are going to study is the Horn-SAT problem Q| where 
all the clauses have at least one negated variable. It is an easy problem and solved in a 
polynomial time. However, notice that as long as K is small fluctuations paly an important 
role in our problem and this could give rise to significantly different behaviors for the 
problem. There is also another problem called 2 +p — SAT j3j which by tuning p goes from 
a 2 — SAT to a 3 — SAT problem. It is close to what we like to do in this paper but here 
we are able to study the easy-hard crossover for general K and it is a more general problem 
to this end. 

In the following we first give the exact solution of 1-SAT problem by a statistical mechanics 
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approach. We find the average number of unsatisfied clauses and the average number 
of solutions in the ground state of the system and explain the origin of their behaviors. 
Utilizing the cavity method and assuming the replica symmetry, we derive a relation for 
the critical point of K-SAT problems. It is found that in general a c oc as p — ► 0. 

Next we obtain the free energy and the distribution of effective fields with the aid of replica 
method and in the replica symmetry approximation. Finally we resort to the numerical 
solution of survey propagation equations for the case of K — 3 and compare the 
extracted critical points with the predictions of replica symmetry assumption. It is found 
that for p > p* ~ 0.17 the replica symmetry breaks at some point ad < a c whereas for 
p < p* we are always in the easy -SAT phase if a < a c . 

The paper is organized as follows. In the next section we define the problem. Section 
IHII is devoted to the study of 1-SAT problem. Assuming the presence of replica symmetry 
we give the results of cavity and replica methods in sections Upland El Survey propagation 
equations for the case K = 3 are numerically studied in section IVI1 Section IVHI includes 
the conclusion remarks of the paper. 

II. THE PROBLEM DEFINITION 

We take iV logical variables {xi\i = 1, . . . , N} where Xi — 1 if the corresponding variable 
is true and otherwise Xj = 0. Alternatively we can speak of N Ising variables Si := 2x» — 1. 
On the other side we have a formula which consist of M clauses which have been joined to 
each other by logical AND. Each clause in turn contains K logical variable selected randomly 
from the list of our N variables. These variables, which join to each other by logical OR, are 
negated with probability p. One obtains the original random K-SAT problem by choosing 
p = 1/2. Here there is an example of a 2- SAT formula with 4 clauses and 10 logical variables 

F := (x 5 V x 7 ) A (x 2 V x 9 ) A {x\ V x±) A (x e V x 10 ). (1) 

A very useful concept in these problems is the factor graph, a bipartite graph of variable 
nodes and function nodes. In FigGJwe have shown the factor graph of formula given above. 
In this figure the logical variables and the clauses have been represented by the circles 
(variable nodes) and the squares (function nodes) respectively. An edge in this graph only 
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FIG. 1: Factor graph of the formula given in EqQ 



connects a variable node to a function node and its style gives the nature of the logical 
variable in the associated clause. Here dashed edges are used to indicate that the negated 
variable enters the clause. We can summarize the factor graph in matrix Cmxn with elements 
C a ,i £ {0, +1, — 1}. In fact C a ^ is +1 or —1 if clause a contains variable i or its negated 
respectively. Otherwise C a ^ = 0. 



III. THE SIMPLE CASE OF K = 1 

We start by g iving the exact behavior of 1-SAT problem by keeping a statistical mechanics 
approach 0, [l(| . We define the energy of a formula as the number of violated clauses, that 
is 

M N 

E[S,C] :=£[l-£C ,^]/2, (2) 

a=l i=l 

where S denotes the configuration of Ising variables. Note that by definition 

M 

Ca,i — U — fi, (3) 

a=l 

where ti and fi give the number of full lines and dashed lines respectively emanating from 
variable node i. The set of {tj, fi} only depends on the structure of the factor graph. 
Utilizing the above facts and summing over spin configurations the partition function reads 

Z[C] := £ e^ E ^ = J] (e~ pti + e~ pfi ) , (4) 

S i 



where (3 = 1/ksT. The free energy per variable, — (1//3JV) ln(Z[C]), still depends on the 
structure of the factor graph and we should take an average over this kind of disorder. 
The probability to have the set {ij, fi} is given by 

P[{tj}) = (m/N M )u(^^y (5) 

Then the averaged free energy in the thermodynamic limit reads 

oo 

/:=7M = a/2-l//9{ln(2) + e- a £ J m {qa)I n {a) m(cosh(/?(m + n)/2))}, (6) 

n,m=— oo 

where we have defined g := 1 — 2p. Moreover J m {a) and J n (ot) are the Bessel functions of 
first kind. Now from Eq.(JH]) one can easily find the average energy per variable 

oo 

e := < E >/N = a/2 - e~ a £ J m {qa)I n {a){m + n)/2 tanh(/3(ra + n)/2). (7) 

m, n=— oo 

In the same way the entropy per variable is given by 

oo 

s := < S >/N = ln(2) + e' a £ J m (qa)I n (a){\n(cosh((3(m + n)/2)) (8) 

m,n=— oo 

-P(m + n) /2 tanh(/5(m + n) /2)}. 

We are interested in the ground state properties of the problem and to this end we need 
to take the limit (3 — > oo in the above relations. After some simplifications we find for the 
ground state energy 

ec = (1 -p)at - e~ a {m + n)J m (qa)I n (a). (9) 

m+n>0 

In Fig |21 we have shown the behavior of this quantity versus a and for some values of p. As 
expected the problem is always in the UNSAT phase where an infinite number of clauses 
are not satisfied in the thermodynamic limit. Indeed the probability that a new clause is 
not satisfied is given by 2p(l — p)a and this quantity has always a nonzero value as long as 
p and a are nonzero. 
The entropy of system at zero temperature is found from Eq.(JBJ) 

SG = e- Q / (a v /l-g 2 )ln(2). (10) 

We have numerically computed this quantity and have shown its behavior with a in FigCJl 
. Clearly the number of solutions is the same for different values of p and small a. It is 
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FIG. 2: Ground state energy of 1-SAT problem from top to bottom for p = 1/2, p = 1/4 and 
p = 1/6. 

due to the fact that for small a each clause constrain a variable regardless of the nature of 
variable in that clause. However for large a the situation is different. In fact when a variable 
contributes in different forms in a number of clauses, with a larger probability it can take 
the values 0, 1 without changing the number of violated clauses. This is the reason for the 
smaller entropy of smaller values of p in FigCJl 

IV. CAVITY METHOD 

In this section we apply the cavity method to our problem. For simplicity we work 
only in the replica symmetric scheme and zero temperature. 

Consider the variable node So with z = 3 neighbors in Fig0] It is known that the cavity 
field experienced by variable So is 

h = ±u(ht\...,h^ ) _ 1 ), (11) 

a=l 

where 

K-l 

u(h u h K -i) = C J] 9(-C r h r ). (12) 

r=l 

Here 6{x) is the known step function which is 1 for x > and otherwise. Note that u and 
h are stochastic quantities. In the replica symmetric framework the probability distribution 





FIG. 4: So experiences a cavity filed determined from the fields experienced by the other neighbors 
of function nodes 1, 2, 3. 



of these quantities are given by [11 1 



Q{u) = cq5{u) + c-5(u + 1) + c + 5{u — 1), 

— 1 oo 

P(h)= P-WSil-^+PoSi^ + Y^P+WSil-h). 



(13) 



l=—oo 



1=1 
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Then it is an easy exercise to show that for h > the probability distribution of h is 



P + (h) = E /*«(*) £ 
z=h r=0 



Z 

2r + h 



2r + h 
r 



c z - 2r - h c r _c r + h . (14) 



Here [. . .] denotes the integer part of enclosed quantity and fica(z) is the degree distribution 
of variable nodes 

f K a(z) = (z Ka /z\) e~ Ka . (15) 

Similarly for h < 0, P-(—h) is obtained by interchanging the role of c_ and c + in the above 
equation. Simplifying the above relations one finds that 

P = e - Ka ^I (2Ka^I^), (16) 

and after some straightforward algebra for P + = Y^hLiP+(h) one obtains 

POO i 

P + + P = e~ Kac + / J Kac+t). (17) 

J Kac- 

Eqs. fllfiJ) and (|T7j) are two independent relations which along with the normalization con- 
dition P + P + P + = 1 determine P(h) and Q(u). However we still need to derive the 
relations between {co, c_, c + } and {Po, P_, -P+}. To this end note that from Eq. (fT2j) one has 

c = l-[pP + + (l-p)P_] K - 1 :=1-Cp, c_=pc p , c + = (l-p)c p . (18) 



Summing the above relations one can use Eqs. (jl6j) and ()17)1 to find the following equation 

c P = [g(c P )} K -\ (19) 



for c p 



where 

S(cp) := 1 - p - (1 - p)P - (1 - 2p)P+. (20) 
For p = 1/2 we recover the known relation for the effective filed distribution 

P(h)=e- Kac "I h (Kac p ), (21) 



where now c p satisfies 



.l- e -«oc Pl ^ K ^ 



2 

Returning to our general problem we find from Eq. (j20j) that g(0) =0 and 



(22) 



dg(c p ) 
dc p 



= 2Kap{l-p). (23) 
8 




FIG. 5: a c versus p for K = 3 and in the replica symmetric approximation. The line shows a power 
law of exponent 2 

Thus as expected Eg. suggests a continues transition for K = 2 and discontinuous 
transitions for K > 3. Indeed for K = 2 the critical value of a is given by 

"c = — 7T~ (24) 

Due to the absence of replica symmetry breaking the above results are exact for K = 2. As 
expected for p = 1/2 we get the known value of a c = 1 and a c —>■ oo for p — ► 0, 1. 
What can be said about a c for general Kl Let us focus on the behavior of a c as p — >• 0. 
First note that just above the critical point P + takes a finite value and from the definition 
of c p , Eq.([T%j). one obtains c p oc On the other hand in g(c p ), c p always appears along 

with a as ac p . Expanding g(c p ) one finds that only for a finite ac p the two sides of Eq. (jl9|) 
would have the same scaling with p. This in turn suggests that the critical value of a should 
scale as p"^^. 

In FigEl we have solved Eq. (|19j) numerically for K = 3. Indeed for K > 3 the replica 
symmetric predictions provide an upper bound for a c For K = 3 we found that as 

expected a c approaches to infinity like p~ 2 as p — > 0. 

V. REPLICA APPROACH 

In the following we will keep the same lines as Ref.^^ to calculate the free energy of 
biased random K-SAT problems in the replica formalism. As before we can write energy or 
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the number of violated clauses as 



m / n \ 

E[S,C]=Y,S + K . (25) 

1=1 \i=i / 



Our goal is to find ln(Z[C]) and to this end we need to obtain 



Z n :=Z\C¥= E e-eZ>: =1 Els:c\, (26) 



The overline denotes averaging with respect to the random structure of the factor graph. 
Due to the independent nature of clauses on can write 

Z n = E C M (27) 

where 



Then for a single clause we can use the fact that 

V i I j=i 

to write the following expression for ( 

C = ^- v E P U[C] (1-P) K ~" [C] E e- pE ^i S ^ +C3 )(l + 0(l/N) + ...), 



(29) 



C , i,...,C Js: =±l h,...,i K =^ 



\ K 1 



(30) 

where z/[C] is the number of minus elements in the set {Cj\j = 1, . . . , i^}. Let us also define 

AT 

i=i 

— * 

where a and Si are vectors of n Ising elements in the replica space. Then apart from some 
irrelevant terms and constants we obtain 

C= ]T p -\P\ {l _ p) K-,[C] x {-C 1 a 1 )...x{-C K a K )e-^^U 8 ^- 1 ). (32) 

C*i,...,Ck=±1 3\,...,3k 



and thus 



Z„~ 5>-"'<*>, /» :=-aHW(x)) + J2x(a)Hx(a)\ (33) 

x(<j) I? 
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where a := M/N and 



K 



W{x);=Y,P v {l-P) K ~ v E n>(^i) II xi-a^e-^^ 5 ^' 1 ) . (34) 

v\K ai,...,a K j=l j=v+l 

In this equation J2u\k means a sum over all the selections of v variables from K ones and 
in the same time it orders these selected variables in the beginning of a i^-member list. 
Now we should find a form for x(a) which minimizes f{x). As in previous studies we 
use the following ansatz in the replica symmetric scheme 



x cr 



^mPHnj^). (35) 

Note that for p ^ 1/2 we do not have the symmetry relation x{—&) = x(a) and so P(m) is 

not an even function. 

Applying the above ansatz we find 

K K 

W(x) = Y,B{v,K-p) / I[dm j P(m j )I[A VtK -. v (m), (36) 

u=0 J j=l a 

where B(i/,K;p) is the binomial distribution 

/ 

B(u,K;p) = 



K 

v 



p\\-p) K -\ (37) 



and 

w»)== e nf 11 ? 5 ) n (i^s) ^ict-j-o, (38) 

Doing the sum over ex one obtains 

A^m) = 1 + if* - 1) ft (^P) fl (^) • (39) 

If we Optimize /(x) with respect to x(a) we find that 

x(a) = Ae^W^M (40) 

where 

:= j^T^)' (41) 
d)x(cr) 

and A is determined from the normalization condition. After some calculations one finds 
the following relation for W'(x) 

K K 

W'(x) = Y,B(v,K;p) / n^ i ^(m i )[z/n^-i^-,(H + (^-^)n^--iH] ) ( 42 ) 

u=0 j=l a + a- 
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where Y\ a± denotes a product over the indices a for them a a = ±1. Moreover in this equation 
K-^Am) = 1 + (e-> - 1) n C-^p) ff Oqp) , (43) 

3=1 3=f 

" / 1 + rrin \ (\ - m. 



^_,_ 1 (m) = i+(e-^-i)n pp n 

3=1 V 7 3=^+1 

We are interested to limit n — > where a; (a) can be written as 

x{a) = J dmP{m)e uln ^\ (44) 

Now doing the standard algebra we find the following self consistency relation for P(m) 

P ( m ) = r due -^m)e- aK+ " w '(™\ (45) 

1 — m 2 J-oo 

where 

K-l K-X 

W\u) = B{v,K;p) / J] dm j P{m ] )[ve uXri y A '- 1 ' K ' Am) ) (46) 
Similarly for the free energy we find 

K „ K 



Finally let us consider the limit (3 — > oo of Eqs. 
with effective fields z given by m = tanh(/?^/2) 



(5f = - \n{2) - a{\- K)Y^B{v,K;p) ( f[ d mj P( mj ) In (A V , K -*) (47) 

^=o 7 j=l 

-a/2Y,B{v,K;p) / J] dm 3 P{m ) [v In (A^i,^) 

+ {K-u) In (A„,k-_„_i)] + 1/2 / dmP{m) In (l - m 2 ) . 

Eqs. fl47|) and (|43j) return the known relations for p = 1/2 ^(J when P(m) is an even function. 

and (J47|) . To this end we should work 
Then for (3 = oo we get 

R( z ) = e~ aK I " due- iuz e aW ' (iu) . (48) 

J — oo 

Now the relation that gives W'(u) reads 

poo r0 

W\u) = K-K[p dzR(z) +{l-p) / dzRiz)}*- 1 (49) 

JO J— oo 

A'-l 

£ B(v,K;p)[v j D v ^ K . v e~ u ^(1,,,.,,^-,,,-*,) 
+ (K-z/) /" D v ^ v ^e u min(i >Zlr .., z ,-^ + i > ...,-^_ 1 )]_ 
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i/=0 



where 



Du-i,k-u := / H dzjR{ Zj ) / H dzjRfa), (50) 

roo v rO A '-! 

D V ^ K _ V _ X := \ JJdzjRizj) ]J dzjR(zj). 

Jo j=i J -°°j= v+ x 



For the free energy in this limit we have 

K 



f = a{l - K) Y B(u, K; p) / D U)K „ V min(l, z 1: . . . , z Vl -z v+1 , . . . , -z K ) (51) 

K 

+a/2Y,B(u,K;p)[u / D v - ltK - v min(l, z x , . . . , -z v , . . . , -z K ^i) 
+ (K - u) J D^k^x min(l, z u . . . , z v , -z u +i, • • • , -z K -i)\ 

/0 roo 
dzR(z)z - / dzR(z)z\. 
-oo JO 

Considering the simple case of K = 1 the effective field distribution reads 

oo / oo \ 

R(z) = e- a /oM/ m («) + E(- 1 ) n/ 2nM[ J 2n- m (a) + / 2 „ +m («)] \6{z-m), (52) 

m=— oo \ n=l / 

which for p = 1/2 returns 

oo 

R(z) = e- Q Y, 4(«)<W (53) 

n=— oo 

Compare the above relation with Eq. |^ which gives the effective filed distribution in the 
cavity method and in the replica symmetric approximation. In fact the two distributions 
are the same as they should be as long as we use an ansatz in which the effective fields take 
integer values. 

VI. SURVEY PROPAGATION EQUATIONS 

In this section we study the behavior of 3-SAT problem by means of numerical solution of 
survey propagation equations ^| . Let us first write the general form of these equations. 
We define t] a ^i as the probability that in a state selected randomly from the existing states 
of the problem, the clause a sends a warning to variable i to take the value that satisfies 
it. This warning is sent if the other members of a do not satisfy this clause. We denote by 
V(a) the set of neighbors of a. Then assuming a tree structure for the factor graph we have 

n ™, (54) 

jeV(a)\i 
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FIG. 6: The survey warning rj a ^,i is determined by the set of surveys rjb-^j. 

where the product is over all the neighbors of a excluding i and (j) is the probability that 
variable j does not satisfy clause a. Let us denote by V*(j) the set of clauses that variable 
j appears in them as it appears in clause a, FiglHl The remaining set of clauses are denoted 
by V™(j). With these definitions P" (j) is given by 

m 



where 



n 



n° + IP ' 



(55) 



rr 



n 



[i- n n c 1 

[i- n a-^i)] n a 



n 



n (i 

66^(j)|o 



(56) 



Now starting from an arbitrary configuration for the warnings sent along the edges of the 
factor graph one obtains the new values of r/s from Eqs. (j54l55l56|) and repeat this procedure 
until reaches to a stationary state. It is believed that in the whole region of SAT phase the 
above equations result in the correct solutions of random K-SAT problems Here we 

apply the same procedure to 3-SAT problem to compute S, the complexity of our problems. 
The complexity of a formula is the logarithm of the number of states and reads |l5j 



M N 
a=l i=l 



US,- 



(57) 
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FIG. 7: From top to bottom: the replica symmetry predictions for a c (RS), survey propagation 
predictions of a c (SP) and ad (SP) for K = 3 and N = 10000. The numerical results have been 
obtained for one realization with the convergence limit equal to 0.001. 



where 



s a = io g [ n (nj_* + + n^J - n n?-J 

jeV(a) jeV(a) 

e, = io g [nr + n° + n+] 



(58) 



and 



(59) 



nr=[i- n (i-^)] n o--rh~i), 

aeV-(i) a£V+(i) 

n+ = [i- n (i-^)] n a-^)- 

aeV+(i) a6V-(i) 

n°= n (i-^)- 

o6V(i) 

In these equations V(i) denotes the set of Z{ neighbors of variable node i, V+(i) is the set 
of function nodes in V(i) that have been connected to i by a full line and V-(i) gives the 
complementary subset. 

It is known that E is zero in the replica symmetric and UNSAT phases and nonzero in the 
hard-SAT phase Increasing a one first encounters the replica symmetry breaking point 
at ad where E takes disco ntinuously its maximum value E m . After this stage E decreases 
and finally vanishes at the critical point a c . One can use these properties of E to compute 
ad and a c . 



To solve the survey propagation equations we used the software given in 17J. In FigJ7|we have 
shown the results of this computation for a c and ad and compared a c with the predictions 
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FIG. 8: Maximum complexity of 3-SAT in terms of p. The parameters are the same as FigJTJ 

of replica symmetric case. As the figure shows the behavior of a c with p is qualitatively 
similar to the one obtained with the replica symmetry assumption. The represented data 
have been restricted to relatively large values of p. It is due to the fact that for smaller 
values of p the complexity vanishes and we are not able to identify a c by looking at E. 
In FigJSJwe also showed the behavior of S m versus p. It is seen that around p* = 0.17 the 
maximum complexity vanishes discontinuously. Then it can be concluded that for p < p* 
we have a simple problem as in the regime of easy-SAT phase. 

VII. CONCLUSION 

In summery we studied biased random K-SAT problems in which a variable is negated 
with probability p. This definition enables us to go continuously from easy random K-SAT 
problems to the hard ones. Certainly this can help us in a better understanding of the typical 
complexity of random K-SAT problems. In this paper we gave the exact solution of 1-SAT 
case and the full picture of general K-SAT problems in the replica symmetry approximation. 
From these results, which are exact for K = 2, one can obtain an upper bound for the critical 
value of a c (p, K). We found that a c (p, K) has a power law behavior p~ TK for p — > where 
tk — K — 1. We studied 3-SAT problems with the help of numerical solution of the survey 
propagation equations and found no replica symmetry breaking transition for p < p* ~ 0.17. 
However in contrast to the tricritical point of 2+p-SAT problem we found that in both sides 
of p* the SAT-UNSAT transition is discontinuous. This phenomenon dose not support the 
current belief that hardness of a problem may stem from the discontinuous nature of its 
transition. Certainly it still demands more studies to have a clear picture of the origins of 
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typical complexity in these problems. 



Acknowledgments 

We thank M. Mezard for helpful suggestions in writing this paper. A. R. is grateful to 
V. Karimipour for the introduction of this beautiful field to him. 



Y. Fu and P. W. Anderson, J. Phys. A 19, 1605 (1986). 
S. Kirkpatrick and B. Selman, Science 264, 1297 (1994). 

R. Monasson, R. Zecchina, S. Kirkpatrick, B. Selman and L. Troyansky, Nature (London) 
400, 133 (1999). 

M. Mezard, Science 301, 1685 (2003). 

C. H. Papadimitriou, Computational Complexity, Addison- Wesley, 1994. 
M. Mezard, |co~nd-mat/0212448| 



G. Parisi, cs.CC/0312011 



O.C. Martin, R. Monasson, R. Zecchina, Theoretical Computer Science 265 (1-2), 3 (2001). 
R. Monasson and R. Zecchina, Phys. Rev. Lett. 76, 3881 (1996). 
R. Monasson and R. Zecchina, Phys. Rev. E 56, 1357 (1997). 
M. Mezard and R. Zecchina, Phys. Rev. E 66 (2002) 056126. 

S. Cook, Proc. 3rd Ann. ACM Symp. on Theory of Computing, Assoc. Comput. Mach., New 
York, 1971, p. 151. 

T. Walsh, APES Report: APES-37-2002. 

A. Braunstein, M. Mezard, M. Weigt and R. Zecchina, cond-ma t/0212451 
A. Braunstein, M. Mezard and R. Zecchina, |cs.CC/0 212002, 
S. Franz and M. Leone, J. Stat. Phys. Ill, 535 (2003). 
http:/ /www.ictp.trieste.it~zecchina/SP 



17 



